How do spins interact coherently with ultrashort laser pulses ? 



Helene Vonesch and Jean- Yves BigotQ 
Institut de Physique et Chimie des Materiaux de Strasbourg, 
UMR 7504 CNRS, Umversite de Strasbourg BP. 43, 
23 rue du Loess, 67034 Strasbourg Cedex, France 

We investigate the interaction of femtosecond laser pulses with spins including relativistic cor- 
rections. The time-ordered magneto-optical signals corresponding to a pump-probe configuration 
are calculated in the case of one electron submitted to a magnetic field and evolving in eight levels 
of the fine structure of a hydrogen-like atom. Our simulations explain the origin of the coherent 
magneto-optical response and ultrafast spin dynamics in ferromagnets excited by intense laser pulses 
as recently reported in Ni and CoPt3 ferromagnetic thin films. Our detailed analysis allows iden- 
tifying the respective roles of the coherent spin-photon interaction and spin dynamics unraveling 
recent controversies about the laser induced ultrafast magnetization dynamics. 



Ferromagnetic thin films can be modified on a fem- 
tosecond time-scale by ultrashort laser pulses [THE • This 
observed phenomenon, as well as the non-thermal optical 
control of magnetic order with light pulses [B , is impor- 
tant as it promises major applications in the fields of 
data-storage and time resolved magnetic imaging for ex- 
ample [7 . Understanding the underlying physical mech- 
anisms requires a formal description of the laser in- 
duced ultrafast magnetization dynamics which involves 
the multiple possible interactions between the photons, 
electrons, the spins and the lattice together with the com- 
plexity of the electronic band structures of the magnetic 
materials. 

Several models have been proposed to describe the de- 
magnetization process occurring in the sub-picosecond 
time scale. The first one involves three interacting baths 
at different temperatures corresponding to the charges, 
the spins and the lattice which are out of equilibrium 
[8]. A quantum model, including the effects of exchange 
and spin-orbit interactions, has accounted for spin flips 
in terms of their dephasing after they redistribute in the 
excited states [5rfTT]. An extension of this model includ- 
ing the combined interaction of the laser and spin-orbit 
successfully described elementary spin-flips or so-called 
lambda processes [12]. Modeling the spins dynamics at 
later times, a few picoseconds after laser excitation, has 
reached a consensus. In that case, the spin-phonon inter- 
action prevails and is responsible for the damping of the 
precession of the magnetization in ferromagnets excep- 
tion made when approaching the Curie temperature for 
which long range fluctuations maintain a non-equilibrium 
spin bath for a long time [T2J [H] ■ A revival of the debate 
regarding the origin of the ultrafast demagnetization oc- 
curred with the assumption that spin flips could be due 
to a mechanism similar to the Elliott- Yafet scattering 
of conduction electrons by magnetic impurities [HI [16] . 
Some experiments do not support this model for exam- 
ple regarding the effect of magnetic impurities on the spin 
dynamics in doped ferromagnetic transition metals [17] . 

More recently a new controversy started regarding the 
interpretation of the time resolved magneto-optical re- 



sponse as a signature of the magnetization dynamics 
[T51 [H)] . The main origin of this controversy lies in the 
distinctions that one ought to make between the time de- 
pendent response function (magneto-optical signal) and 
the system's dynamics (the spins populations). In the 
present work we show that there is a straightforward 
manner to clarify the debate by considering separately 
the coherent and population dynamics in the magneto- 
optical response. Indeed it is known that in metals coher- 
ent magnetism is important as shown in a recent detailed 
study of the charges and spins dynamics performed in Ni 
and CoPt3 ferromagnetic films [2"U] . 

Let us consider the simplest possible system consti- 
tuted of eight quantum levels interacting with the laser 
field such that the interaction takes into account the 
relativistic corrections to the quantum electron dynam- 
ics, including the spin-orbit interaction. We determine 
the response function from the density matrix formal- 
ism including the time ordered third order nonlinear 
terms. This approach allows understanding the main dif- 
ferences between the coherences and spins populations. 
We consider a one-electron Hamiltonian with an effective 
Coulomb interaction perturbed by the spin-orbit and ki- 
netic momentum-laser vector potential interactions. Us- 
ing such Hamiltonian, applied to the band structure of 
metals, P. N. Argyres has shown that static magneto- 
optical Kerr and Faraday effects in ferromagnetic mate- 
rials mainly come from spin-orbit interaction with the 
ionic field [21 . For the purpose of understanding the 
main steps of our approach, which primarily aims at de- 
termining the respective roles played by the coherent ver- 
sus population dynamics in the response function (Fara- 
day rotation), we consider the case of a simple discrete 
eight level system representing a reduced hydrogen-like 
system. 

The relativistic contributions to the ultra-fast magneto- 
optical dynamics are considered via the Foldy- 
Wouthuysen transformation of the Dirac equation which 
enlightens the various interaction terms between spins 
and femtosecond laser fields. We add to Argyres' ap- 
proach the terms of the Foldy-Wouthuysen transforma- 
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tion to second order in — for the electron from an 
Hydrogen-like atom submitted to different static fields 
: the ionic field Ej and its associated central ionic po- 
tential Vi(r) and a strong static homogeneous magnetic 
field %e z with potential vector Am = — |RABj\/- The 
electron is interacting with a laser field described in the 
Coulomb gauge associated to an homogeneous electric 
field El, to its colinear vector potential Ai and related 
magnetic field B^. 

The Hamiltonian H corresponding to no interaction 
with the laser reflects the Zeeman splitting in a strong 
magnetic field. At this point the degeneracy of the spin 
states are already lifted and the effect of spin-orbit inter- 
action is to slightly shift some energy levels. 
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Neglecting the second order terms in the laser vector po- 
tential the interaction Hamiltonian is : 
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With the kinetic momentum operator 11= S [R ; H ] : 



n 



qA 



2mc 2 



S A E, 



(3) 



We apply the electric dipolar approximation to the in- 
teraction Hamiltonian and consider the 2s and 3p levels 
of the Hydrogen atom ; their energy difference without 
spin-orbit coupling and static magnetic field is associ- 
ated to the frequency ujq as shown in figure [l] This 
energy is close to the typical photon energies used in 
experiments. We do not consider the magnetic dipo- 
lar interaction term which is off-resonant. The two first 
terms of equation 2 imply a variation of the projection 
of the orbital momentum along the quantization axes 
Al = ±1,AZ 2 = ±l,As 2 = and the terms contain- 
ing the spin (H a and Hp) and spin-flip transitions with 
A; = ±1, Al z = 0, As z = ±1. Table [I] shows the orders of 
magnitude of each interaction matrix element compared 
to the predominant interaction term < j|^p.Ai|i > for 
two states \i > and < j\ among the eight states sketched 
at figure 1. In the case of our hydrogen-like model, we will 
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see that the contributions H a = ^j[SAEj.A £ and 
Hp = 9m 9 2 c 2 [p A S].Ei can be neglected in the magneto- 
optical response but not in the spin dynamics. It is not 
the purpose to discuss here their influence in the case of a 
ferromagnetic solid, where additional phenomena such as 
a dynamical anisotropy can be induced by the presence 
of an external field [22] , the laser field for example in our 
case. 
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FIG. 1: Considered transitions in the Hydrogen-like atom. 
For each level we indicate l z ,s z . a± stands for a circularly 
polarized field along e x ± ie y . 



TABLE I: Scale of the considered interaction matrix elements 
in Hydrogen compared to the canonical momentum-laser vec- 
tor potential interaction matrix element. 
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In order to model magneto-optical pump-probe exper- 
iments, the evolution of the system is calculated in the 
Liouville formalism to the third order of the laser per- 
turbation [23] . The relaxation time of the coherences 
(Pnmjf 1 7^ m ) is T2 whereas it is T\ for the population 
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plicity they are assumed to be the same for each tran- 
sitions. The electric field of the laser pulses is given 



byEi(t) = |[e(t,r)e- 



- e* (t, T)e WLt ]e x where e(t, r) 



is a gaussian centered on the pump-probe delay t or 
in the case of the probe or pump fields. The laser fre- 
quency Wi is equal or close to ujq. The effective dipolar 
moment of the system D gives the first and third or- 
der polarizations of the atom P^'(t) — Tr(p^T)) and 
p( 3 )(i,r) = Tr(p( 3 )D). The polarization is calculated 
using the rotating wave approximation and depends on 
both the time t and delay r between the pump and probe 
pulses. As we do not consider any propagation effect in 
this simple atomistic approach, we set arbitrarily that 
the radiated electric field at order (n) is E^™- 1 = P^ n \ 
We consider pump and probe linearly polarized in the 
plane perpendicular to the quantization axes and cal- 
culate the dynamical magneto-optical rotation induced 
by the sample in the Jones formalism. In order to fol- 
low the common experimental procedure for magneto- 
optical measurements with a polarization bridge, we de- 
fine E' (3) (£, r) as the rotated E (3) (£, r) due to a half- wave 
plate tilted by an angle of ? with respect to the (e^, e^) 
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axes, the measured differential intensities AI x ^ y ) are: 
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The differential rotation for a positive magnetic field 
A8 + b m (t) = [AI X — AI y ] + B M is linked to the first and 
third order field's amplitude and to their instantaneous 
rotations ®+ B (t) and Q^ B (t,r); here we assume the 
ellipticity to be negligible. 



A6 +Bm (r) = «n[0S M (i, r) + 6^ M (t)] (5) 
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The resulting magneto-optical signal is then propor- 
tional to: A9(t) = AQ +Bm (t) - A0_j3 A/ (t ) obtained 
for the two directions ±Bm of the magnetic field. As 
done experimentally, we perform a differentiation on the 
magnetic field and normalize the resulting rotation with 
the linear magneto-optical rotation Q±b m : 
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Considering the time ordering of a single probe 
and two pump pulses one can distinguish three terms 
corresponding to different coupling between pump, 
probe and the polarization of the system as de- 
scribed in the case of charge dynamics by Brito- 
Cruz et al. [25]. First, the time ordering of 
the fields e pump (t = 0)e pump (t = 0)e probe (t = r) describes 
the "population dynamics" induced by the pump 
field c pump ,e* ump . It is maximal for positive pump- 
probe delays and relaxes with the population de- 
cay time T\. The second term named "pump- 
polarization coupling" (PPC) is given by the se- 
quence e* ump {t = 0)e probe (t = T)e pump (t = 0) which cor- 
responds to a convolution of the pump pulse with 
the exponential decay of the density matrix's co- 
herences within the Ti time. The third term 
£ P robe(t = r)e* ump (t = 0)e pump (t = 0) is maximum at 
negative pump-probe delays and corresponds to the co- 
herences generated by the probe which couple to the 
pump, sometimes also named "pump-perturbed free in- 
duction decay" (PP-FID). These two latter terms, also 
named "coherent terms" hereafter, are strongly related to 
the dephasing of the polarization as pump and probe can 
only couple if the density matrix's coherences are non- 
zero. Figure [2] shows the dynamical differential magneto- 
optical rotation calculated for each of these three terms 
using the interaction Hamiltonian in equation 2. 

In the numerical simulations, we have considered gaus- 
sian pulses of width 10 fs and a probe field 10 times less 
intense than the one of the pump. The dephasing time 
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FIG. 2: Magneto-optical rotation for the population 
(black), pump-polarization-coupling (red), pump-perturbed 
free-induction decay (blue) and total (green) signals. (1): 
on resonance (Twjl = huio = 1.89eV) and (2): off resonance 
{fwjL = L51eV). 



of the coherences and the lifetime of the differences of 
population are chosen to be 10 fs and 100 fs. The static 
magnetic field is 1 Tesla. The pump and probe laser 
fields are taken linearly polarized along e x . The unper- 
turbed ground-level populations are 0.9 for s z — =p| and 
0.1 for s z = ±5 for a magnetic field ±Bm- These results 
clearly show that the coherent terms are important in 
the magneto-optical signals and that their contributions 
respective to the spins population dynamics depends on 
the laser detuning. 

The next issue that we address now is the difference 
between the pump-probe magneto-optical signal and the 
spin and orbital momentum dynamics. Towards that 
purpose we determine two different quantities: the pro- 
jection of the spin and the orbital momentum operators 
?( 2 )\ „„a /r( 2 ) 



(Sz ) and (Lg ') along the quantification axis e 2 . Note 
that the dynamics is now represented by the trace of the 
operators multiplied by the second order nonlinear terms 
of the density matrix as they are related to the density 
matrix's populations. The spin dynamics are due to the 
H a and Hp spin-flip terms. Figure 3 shows their dynam- 
ics for the three time orderings of the pump and probe 
fields. The population terms (fig. 3a and 3b) are repre- 
sented as a function of time t as they do not depend on 
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FIG. 3: Spin and orbital momentum dynamics for 
the three pump-probe field sequences, (a): Tr[p^S z (t)] 
and (b): Tr[p^ 2 ' L z (t)] for the population term. (c): 
R{J dtTr\p {2) S z {t,r)]} and (d): Q{J dtTr[p (2) L z (t, r)]} for 
the pump-polarization coupling term, (e) and (f): same quan- 
tities as (c) and (d) but for the pump-perturbed free-induction 
decay term, (c) to (f) have the same normalization factor. 



the probe field. The two coherent terms are integrated 
over t and represented as a function of pump-probe delay 
t, fig.3c and 3d for "PPC" and fig.3e and 3f for "PP- 
FID" terms, as they explicitly depend on time ordered 
sequences involving both the pump and probe pulses. 



(si 2 ') and (Lf) are respectively real and imaginary be- 
cause they are generated by the pulses e* and e 



-(2 
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In order to calculate real quantities one should also have 
considered e pum p and e* robe which does not correspond 
to the experimental configuration chosen here. 

In conclusion, we have shown that in an ultra-fast 
magneto-optical experiment the time ordering of the 
pulses has to be taken into account especially in order 
to distinguish the coherent response from the popula- 
tions dynamics. Out of resonance the magneto-optical 
coherent signal increases. More importantly the spin and 
orbital momentum's dynamics, due to the second order 
in perturbation of the density matrix, both participate 
to the coherent magneto-optical response. Extrapolating 
the present results to more complex magnetic systems 
shows that the coherent spin-photon interaction can be 
used to manipulate the magnetization of spin devices at 
the femtosecond time scale. 
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